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In a GaN/AIGaN field-effect terahertz detector, the directional photocurrent is 
mapped in the two-dimensional space of the gate voltage and the drain/source bias. 
It is found that not only the magnitude, but also the polarity, of the photocurrent can 
be tuned. A quasistatic self-mixing model taking into account the localized terahertz 
field provides a quantitative description of the detector characteristics. Strongly lo- 
calized self-mixing is confirmed. It is therefore important to engineer the spatial 
distribution of the terahertz field and its coupling to the field-effect channel on the 
sub-micron scale. 
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Field-effect transistors (FET) based on high-mobility two-dimensional electron gas 
(2DEG) are promising candidates for room-temperature high sensitivity terahertz detec- 
tors." The related Dyakonov-Shur's non-quasistatic theory for the plasmon-induced photo 
volt age /current has been studied extensively in the resonant and non-resonant regimes.- 
In the former regime, the channel length L is shorter than the plasma-wave decay length 
(wr ^> 1, sr ^> L), where uj is the plasma frequency, r is the momentum relaxation time 
of travelling electrons, s is the plasma wave velocity. A resonant detection requires a short 
electron channel and a special drain/source boundary condition to enhance the Q factor 
of the plasma-wave cavity. The experimentally observed resonant detection is usually re- 
alized in nanometer-sized FETs at cryogenic temperatures.- - - In contrast, a non-resonant 
detection (ur < 1, sr <C L) relies on the rectification/self- mixing of damped plasma waves. 
As has been pointed out in Ref. Ill and |5|, a self-mixing detector can deliver a high respon- 
sivity in a broad frequency range far above the transistor's cut-off frequency, even in a 
long-channel FET at room temperature. In the self-mixing model, a terahertz excitation is 
usually assumed at the edges of the gate with the source being grounded. Ref. |5| studied 
the spatio-temporal carrier dynamics in the electron channel. The simulation showed that 
the plasma wave decays rapidly from the excitation point toward the other end of the 
channel with a decay length less than 50 nm for silicon at room temperature. It is within 
such a short distance where the terahertz photovoltage is built up. The effect of spatial 
distribution of electron density tuned from the linear regime to the saturation regime by 
the drain/source bias (Vds) and the gate voltage (Vq) has been studied both theoretically 
and experimentally-^ However, a complete map of the terahertz response in the 2D space 
of Vds an d Vq is not yet available. It is of great utility to further reveal the spatial dis- 
tribution of the self-mixing, so that a more accurate device model may be constructed for 
detector optimization. In this letter, we report a full characterization of a self-mixer by 
mapping out the photocurrent in a 2D space of Vds an d Vq, covering the linear (LR) and 
saturated regimes (SR) and the transition regime (TR) in between. A technique to reveal 
the inhomogeneous distribution of terahertz field and the resulting self-mixing is introduced. 
A quasistatic self-mixing model is constructed and examined, showing a good quantitative 
fit to the experiment results. 

To focus on the device physics, we present the experiment data obtained in a GaN/AlGaN 
field-effect terahertz detector similar to that reported in our previous work.- We have to 



mention that similar results were also observed in other detectors with different antennas 
or a different channel length.- A scanning-electron micrograph of the detector and the mea- 
surement setup are shown in Fig. QJa) The key features of the detector are (1) a 2-/^m 
field-effect gate is connected to an antenna block (g-antenna), (2) two isolated antenna 
blocks (z-antenna) sit by the gate with an equal gap of 1.5 /im, (3) the drain and source 
ohmic contacts are set about 64 /im away from the gate. At 77 K, the electron mobility 
and the density are /i = 1.58 x 10 4 cm 2 /Vs and n s = 1.06 x 10 13 cm -2 , respectively. The 
width of the electron channel (mesa) is W = 4 /im. More details of the detector have been 
presented in Ref. |8j. The terahertz wave, generated from a backward- wave oscillator (BWO: 
/ = uj/2tt = 903 GHz), is collimated and focused onto the detector cooled at 77 K. The 
polarization of the terahertz wave is parallel to the i-antennas. As shown in Fig. Ha), we 
used a standard lock-in technique, in which the terahertz wave is modulated by a chopper 
at frequency J'm = 317 Hz. The directional photocurrent is calculated from the magni- 
tude and the phase of the AC voltage across a sampling resistor (10 Q). The DC voltage 
source (Vs or Vd) and the sampling resistor are connected to the drain electrode and the 
source electrode, respectively, or vice versa. The order of the lock-in amplifier's (LIA-S 
or LIA-D) differential inputs (A and B) is kept consistent with the direction of a positive 
photocurrent. The detector is first characterized without terahertz irradiation. As shown in 
Fig. QTb), the current- volt age curves exhibit the standard characteristics of a FET. Accord- 
ing to the gradual-channel approximation (GCA), the channel current in the linear regime 
(LR: Vds < Vq — Vt) is written as 



e^Wn x — , (1) 



where, the channel potential V x varies from to Vos from x = to x = L, and Vp is the 
threshold gate voltage. The local electron density is controlled by the effective local gate 
voltage: n x = CqVgeI^ = Cg(Vg — Vt — V x )/e, where Cq is the effective gate-channel 
capacitance. By measuring the differential conductance (Go = &Ids/&Vds ~ ^WCg(Vg — 
Vr)/L) at Vds = V as shown in Fig. H](c), the threshold voltage and the gate capacitance 
are estimated to be Vt ~ —4.34 V and Cq ~ 0.4 /xF/cm 2 , respectively. In the GCA model, 
the drain-source current and the local channel potential can be written as 
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fiC G W [2{V G - V T )V DS - VI s ] /2L if V DS <V G - V T , 

^C G W [(V G - V T ) 2 + A( V DS -V G + V T )] /2L if V DS > V G - V T , 



(V G - V T )[1 - (1 - x/L LR ) 1 l\x= [0, L] if V DS <V G -V T , 
(V G - V T )[1 - (1 - x/L SR y/%x= [0,L S r] if V DS >V G - V T , 
where, parameter A ~ 0.08 V describes the degree of the effective channel-length modulation 
(L -> L SR = L/[l + X(V DS -V G + V T )/(V G - V T ) 2 }) in regime SR, and L LR = L(V G - 
Vt) 2 /[2(V g — Vt)Vds — Vds} i n regime LR. From Eq. [2] and the experimental data shown 
in Fig. QJc), the charge density and its derivative dn/dV G E at different locations can be 
numerically calculated. 

Upon terahertz irradiation with a frequency of u and an energy flux of Pq, both a hori- 
zontal (£, X E ) and a perpendicular (C, Z E ) terahertz field in the channel are induced. Here, 
the free-space field strength Eq is determined from P Q and the free-space impedance Z by 
Po = Eq/2Z . Parameter £ x and £ z stand for the field enhancement factors for the horizontal 
and perpendicular field, respectively. As shown in Fig. Et^a,b) from a FDTD simulation, the 
horizontal field is concentrated in the gaps between the gate and the i-antennas, while the 
perpendicular field is mainly distributed under the gate. Both fields are stronger at the 
source side of the gate than at the drain side. The horizontal field vanishes and changes its 
phase by 180° at x c ~ +0.3 /xm, while the perpendicular field keeps its phase constant along 
the channel (Fig. [2](c,d)). Since the terahertz fields are highly localized within a distance of 
200 nm to the edges of the gated area and the phases evolve in a continuous manner, the 
effective electron channel can be modeled as a chain of nanometer-sized FETs excited by a 
spatially decaying terahertz field. At 77 K, the momentum relaxation time is about r ~ 2 ps, 
which corresponds to a mean-free-path of < 1 /im. In spite of the fact that ut ps 10, in 
each nanometer-sized transistor, the electron transit time (td <C 200 nm/% ~ 2 ps) and 
the plasma-wave transit time (r p <C 200 nm/s w 0.2 ps) become smaller than the period 
of the oscillating field, i.e., ojto/ p <C 10. In the following, we will present the map of ter- 
ahertz response in the 2D space of Vds an d V G and verify a quasistatic self-mixing model 
taking into account the non-uniform terahertz excitation. Under the quasistatic assump- 
tion, the local effective channel potential and the effective gate voltage can be written as 
V x — > V x + £, x E cos(ujt + 4> x ) and V G — > V G + £ Z E cos(ut + <p z ) with d^/dx = £ x and 



d^ z /dz = £ z , respectively. By substituting the above quantities into Eq. [Hand integrating 
along the channel, we obtain the DC terahertz photocurrent 
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where <fi = <p x since <p z = 0. Current i xz ex £ x !; z cos0 is induced by both horizontal and per- 
pendicular fields, while i xx oc £ x £ x is only from the horizontal field. For numerical simulation, 
we approximate the term i xz as 

-Z P z / ■j=—£ x £ z co8<l>dx, (3) 



2L J dV GE 

where, z = £ z /£ is the effective distance between the channel and the gate. Due to the phase 
change at x c , mixing term i xz changes from a positive value at the left side to a negative 
value at the right side of x c . In the case of the antenna design shown in Fig. [TJ term i xz is 
the major contribution to the induced photocurrent.— 

In Fig. [3(a), the photocurrent as a function of Vds and Vq is mapped in a color-scale 
plot. For comparison, a simulated map based on the above quasistatic model is shown 
in Fig. [3](b). In the simulation, we used two sets of input data. One of them includes 
the parameters such as the electron mobility, the gate geometry, the threshold voltage and 
the gate-controlled electron density obtained from Fig. QJc). The other set of input data 
includes the spatial distribution of the terahertz fields obtained from the FDTD simulations, 
as shown in Fig. [2j The simulation reproduces the main features of the polarity change and 
the magnitude variation. In regime PR, the whole channel is pinched off and the terahertz 
mixing is disabled. In the linear regime, there is only a very weak photocurrent. A strong 
photocurrent is produced in regime TR and SR, and strongly depends on the applied bias 
and gate voltage. When the bias is applied at the drain, a negative photocurrent is flipped 
into a stronger positive photocurrent when the gate voltage is swept from regime TR into 
regime PR. In contrast, when the bias is applied at the source and the gate voltage is 
swept, a positive current is flipped into a weaker negative current. When the bias is swept 
and switched to either the drain or the source with a constant gate voltage at -4.1 V, the 
photocurrent differs in polarity and its magnitude becomes fairly independent of the bias 
voltage. 



The spatial distributions of the terahertz fields (£ x £ 2 cos <j)) and the derivative of the 
charge density (dn/dV GE ) at V G = —4.1 V and V D = 2.0 V or V s = 2.0 V are plotted in 
Fig. E(a). Due to the diminishing of charge modulation at the side where a positive bias is 
applied, the local self-mixing is suppressed. Hence, a polarity change is observed between 
the cases of Vd = 2 V and Vs = 2 V. To further illustrate the terahertz responses in different 
regimes, we mark three lines (Vs — V, Vd = Vs — V, Vd — 2 V) on the maps shown 
in Fig. [3j Along the line at Vd = 2 V, we select five conditions labeled as a, (3, 7, 5, and 
X, corresponding to Vq = —1.0, —2.0, —3.0, —3.5, and -4.1 V, respectively. The location- 
dependent charge modulation is plotted in Fig. HJ^b). By decreasing the gate voltage from 
a to x 5 the charge modulation becomes stronger and the location of the maximum shifts 
from the drain side to the source side. This evolution gives rise to the observed transition 
from a weak negative current to a strong positive current along the line at Vd = 2.0 V. In 
Fig. Etc-e), the experimental gate-controlled photocurrent is compared to the simulation. 
At zero bias, the photocurrent resembles what we observed in our previous studies.- 1 ^ At 
Vd = 2.0 V, the simulation fits well to the experimental data. As shown in Fig. 11(e), when 
the bias is applied at the source side, a transition from positive current at Vg ~ — 2 V 
to a weaker negative current at Vg ~ —4.1 V is reproduced. There is, however, a strong 
suppression of photocurrent from Vg = —3.5 V to Vg = —2.5 V (also marked as vertical 
double arrows in Fig. [3]), while according to the quasistatic model a strong positive current 
should occur. The absence of photocurrent in region SR is also observed in our detectors 
and other detectors with symmetric antennas.— The common feature is 
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that the self-mixing in regime SR is suppressed as long as the electron density at the same 
side of the strong terahertz field is set to be mostly depleted. The deviation between the 
experiment and the model indicates that there is a different transport mechanism. The most 
likely candidate is the plasma- wave induced photocurrent. To reveal the physics behind, we 
will conduct systematic experiments at lower temperatures on different antenna designs to 
tune both the amplitude and the phase of the terahertz field. Accordingly, a non-quasistatic 
model taking into account the localized terahertz excitation will be examined. 

In conclusion, the terahertz response in a GaN/AlGaN field-effect detector was mapped. 
We directly probed the local self-mixing by selectively depleting the 2DEG near the edges 
of the channel. A quasistatic model was constructed and well describes the main features 
of the terahertz response. The finding confirms that non-uniform terahertz fields and hence 
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self-mixing are induced in the channel. The model provides valuable information as to 
how to construct a high responsivity terahertz field-effect detector. The deviation between 
the experiment and the model suggests that a full model taking into account the dynamic 
transport and the spatial distribution of terahertz excitation is desirable. 
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FIG. 1. (a) Measurement circuit diagram. The inset is an artificially-colored scanning-electron 
micrograph of the detector, (b) The I-V characteristics of the detector, (c) The differential 
conductance at Vos = V as a function of Vq- (d) A color-scale plot of the differential conductance 
as a function of Vq and drain/source bias (Vb or Vs) which is applied either at the drain or at 
the source. The dashed lines separate the map into the linear regime (LR), the transition regime 
(TR), the saturation regime (SR), and the pinch-off regime (PR). 

FIG. 2. Spatial distributions of the field enhancement factors (a, b) and phases (c, d) of the 
horizontal and the perpendicular field. 

FIG. 3. The measured (a) and the simulated (b) terahertz response as a function of Vds an d Vq- 
Three vertical lines are marked at Vs = 2 V, Vd = V5 = V, and Vd = 2 V. Along the line at 
Vd = 2, five bias conditions are marked by a, j3, 7, 5, and x, corresponding to Vq = —1.0, —2.0, 
—3.0, —3.5, and —4.1 V, respectively. In the SR regime when the bias is applied at the source, 
there is a clear deviation between the experiment and the simulation. 

FIG. 4. (a) The spatial distributions of the charge modulation dn/dVcE and the self-mixing term 
£,x£,zCos(j) at Vq = —4.1 V with a DC bias of 2 V applied at the source or at the drain, (b) The 
spatial distributions of charge modulation at different gate voltages and at a constant drain bias 
(V D = 2.0 V). The mixing current as a function of V G with (c) V D = V s = V, (d) V D = 2.0 V, 
and (e) V s = 2.0 V. 
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